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Abstract 

We present measurements of the pion transverse momentum {pt) spectra in 
central Si-nucleus collisions in the rapidity range 2.0 < y < 5.0 for pt down 



to and including pt = 0. The data exhibit an enhanced pion yield at low 

Pt compared to what is expected for a purely thermal spectral shape. This 

enhancement is used to determine the A resonance abundance at freeze-out. 

The results are consistent with a direct measurement of the A resonance yield 

by reconstruction of proton-pion pairs and imply a temperature of the system 

at freeze-out close to 140 MeV. 
PACS number: 25.75. -Fr 
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Collisions of heavy nuclei at ultrarelativistic energies produce a zone of hot, compressed 
matter. Information from measurements of transverse energy production |jl|,|^ and baryon 
distributions IQ-Q indicate that baryon densities up to ten times normal nuclear matter 
density are reached during the collision This highly compressed system then expands 

1^ until its constituents cease to interact, i.e. "freeze out". The expansion is reflected in 
the slopes of transverse momentum spectra at midrapidity, which systematically become 
flatter with increasing particle mass |To|JTT[| . Very recently [T^ sidewards flow was directly 
identifled for Au+Au collisions at AGS energy. Such flow effects demonstrate the presence 
of large thermal pressures, and should provide information on the equation of state of the 
hot and dense matter formed in the collision. At the same time, the connection between the 
transverse momentum spectra of hadrons and the temperature of the flreball at freeze-out 
is complicated by additional parameters such as flow velocities and flow proflle. 

To provide information on the composition of the flreball formed in the collision, and to 
get an independent measurement of the freeze-out temperature we report here measurements 
of the double differential cross sections for charged pions near the beam rapidity {y = 
3.4) in central ^^Si+Al, Pb collisions at piab = 14.6 GeV/c per nucleon. The aim of these 
experiments was to measure transverse momentum (pt) spectra for pions in a kinematical 
region where predicted enhancements at low pt would show up strongly. The experimental 
status of low Pt phenomena in high energy nuclear collisions and possible interpretations 
have been summarized recently [|1^. At AGS energies a major source of the enhancement is 
expected to be [ll4| , p!5| the pions produced by the decay of the A(1232) resonance. The decay 
preferentially populates the spectrum at low pt wherefrom the abundance of the A at freeze- 



out can be inferred. This can then be used [T^ to determine the system's true temperature. 



To provide further support for the feeding scenario and independent information on the A 
abundance, we also present the results of an effort to directly reconstruct the A++ using the 
pvT^ invariant mass spectrum. 

The experiment was performed using the E814 apparatus at the AGS at Brookhaven 
National Laboratory. The apparatus is described in detail elsewhere [ITH^,]!?!]. A 14.6 GeV/c 



per nucleoli ^^Si beam was incident upon Pb targets of thicknesses of 1.1 and 2.2 g/cm^ and 
Al targets of 0.33 and 0.66 g/cm^, corresponding to 1.2 and 2.4 % of a silicon interaction 
length, respectively. Collision centrality was determined via a charged particle multiplic- 
ity measurement in the interval 0.85 < r] < 3.8. Experimental details and the connection 
between centrality and charged particle multiplicity are discussed in []TB| . 



Particles emitted in the forward direction were accepted and analyzed by a forward 
spectrometer. We define z along the incident beam, y vertically upward and x so as to 
make a right handed coordinate system. The spectrometer aperture — 115mr < 6^ < 14mr 
and l^^yl < 21mr was defined by a Pb/steel collimator. Accepted particles pass through a 
dipole magnet and are momentum analyzed via a pair of drift / pad tracking chambers. The 
momentum resolution of the spectrometer has been modelled using the GEANT package 
|[T9| . At the low field used in the present measurements, the resolution in momentum p is 
dominated by multiple scattering and is nearly uniform in momentum with 6p/p ~ 4.1% 
for the momentum range considered here. Scattering in the target creates a distortion in 
Pt without significantly altering p. This effect together with all other imperfections implies 
that 6pt < 4 MeV/c for pt < 100 MeV/c and Spt/pt = 4% for larger pt values. All data are 
presented in 10 MeV/c pt bins. 

Time-of-Flight and hence velocity is determined by one of two scintillator hodoscopes 
located 12 m (200 psec resolution) and 31 m (350 psec resolution) from the target. The 
spectrometer is capable of separating protons and pions up to p = 7 GeV/c. Background in 
the pion sample due to kaons and unrecognized decays is less than 10%. Possible electron 
contamination was investigated in two ways. First, for momenta below 0.5 GeV/c electrons 
can be separated via time-of-flight. In this momentum range, the electron to pion ratio is 
observed to decrease with increasing p and is close to 5% at 0.5 GeV/c. Simulations imply 
that the primary source of electron contamination is photon conversion in the target. A 
comparison of results for the 1% and 2% targets shows no statistically significant evidence 
for electron contamination in the pion sample presented here. We conclude that electron 
contamination can be neglected. Consequently, the data obtained with the 1.1 and 2.2 % 



targets were combined. 

Fig. |1| shows a summary of the measured vr" transverse mass (m^ = \Jpi + m"^) spectra 
for central {u /cfgeo < 2%) Si+Pb colhsions. Here m^r is the pion rest mass. The vertical axis 
is X d"^ N / drritdy , the representation in which a Boltzmann (or thermal) distribution 

is a pure exponential in rrit. Our acceptance in rrit is largest for rapidities 3.0 < y < 4.0. 
The measurable range in rrit is limited by the geometrical opening of the spectrometer at 
low y and particle identification at high y. The solid curves in Fig. show Boltzmann 
fits to the data, with the fit interval restricted to {rrit - f^i-w ) > 160 MeV/c^. The data 
exhibit a significant enhancement over this functional form at low rrit. This is most clearly 
seen in the inset which shows the ratio of the data to the thermal fit on a linear scale. 
The dashed line demonstrates that the data rise also faster than expected from an invariant 
distribution l/rrit x d'^N/drritdy oc exp(— m^/T). Results for 7r+, and for the Al target, (not 
shown here) exhibit enhancements of similar strength as those in Fig. ||. The strength of 
the enhancement does not vary significantly with centrality in the range cr/ageo = (2-10)%. 

The large stopping at AGS energies implies that the fireball formed in the collision is 
baryon-rich. Therefore, baryonic resonances and in particular the A(1232), in the following 
simply denoted A, with its low decay momentum of 227 MeV/c are therefore anticipated to 
be the major contributors to pion spectra at low pt [|l^,|l^. One can extract from the pion 
spectra an estimate of the fraction of pions resulting from A decay. To do this we have com- 
puted the pion spectral shape by superposition of direct thermal pions (with spectral slope 
adjusted at high rrit) and pions from A decay for various ratios / = n a /tt direct- Rapidity 
and Pt distributions of the A were taken to follow predictions by the RQMD model pO[ but 



the assumption that the A distributions follow those of the measured protons yields very 
similar results. Fig. shows the ratio of the pion spectra to the fitted Boltzmann distribu- 
tion. The data are compared to model calculations using /=0.4 and 0.6. The calculations 
bracket the data quite well and establish that the ratio / = 0.5 ±0.1 implying that 1/3 of 
all pions come from A decay without further interaction. The ratio between the number of 



pions and nucleons for central Si+Pb collisions is observed to be {7T/N)exp ~ 1-07 ITTI , ^ 



Hence the fraction of nucleons excited to the A resonance at freezeout is 0.36 ± 0.05. 

The measured pion spectra are also well reproduced by cascade models such as RQMD 
pO| and ARC |^ which have the A resonance explicitely built into the collision dynamics. 
This is illustrated by the dashed lines in Fig. |I] representing the RQMD prediction, which 
accounts for shape and absolute yield of the data once the experimental trigger conditions 
are incorporated. Predictions using the RQMD and ARC models are in good agreement 
with all our pion data (for both charges and both targets). The dominant source of the rise 
at low rrit in pion spectra calculated with RQMD can be traced back [^] to the A(1232) 
resonance decay. The overall predicted freeze-out A excitation probability of 0.35 is very 
close to the experimental value given above. 

Close inspection of Fig. |l] reveals that in addition to the low pt enhancement discussed 
above there is also visible in the data an increase with even steeper slope at very low 
transverse momenta {pt < 50 MeV/c corresponding to rrit — m.^ < 0.01 GeV/c^). A similar 
effect is seen |^ in our recent measurement covering backward rapidities with a different 
detector. Whether this is also due to resonance decays (such as t] decay, see [|T3|) or has a 
more exotic origin (such as chiral symmetry restoration [^) remains to be quantitatively 
explored. 

To get an independent measurement of the number of A's at freeze-out in Si+Pb collisions 
we have reconstructed the A++ via its decay to p7r+. This is the most easily measured 
of all A decays since (i) its branching ratio is nearly 100%, (ii) all particles in the final 
state are charged, and (iii) there is no interference from A decay (which could disturb A*^ 
measurements). Additionally, the asymmetry of the E814 spectrometer makes it best suited 
for like-sign pair measurements. 

The invariant mass for pvT^ pairs was reconstructed for central {cr/ageo < 10%) collisions. 
Protons near beam rapidity (y > 3.1) were rejected from the sample since they are in part 
projectile fragments. Pions below rapidity 3.0 were rejected since the A decay kinematics 
does not permit such pions into the E814 spectrometer acceptance with the proton rapidity 
cut used. 
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Fig. 1^ shows a summary of our measurement of the A^"*" using pvT"*" pairs. The analysis 
employs the "mixed events" technique. This method allows to separate a small signal from 
a large combinatorial due to uncorrelated pairs. One determines the shape of the combina- 
torial background by constructing an invariant mass spectrum using (uncorrelated) protons 
and pions from different events. The resulting distribution is normalized (see below) to the 
true pair spectrum and subtraction yields the signal. To normalize, we treat the combina- 
torial spectrum as a function of a single free parameter, the normalization constant. This 
function is fitted to the high mass end {Minv > Mi) of the true pair spectrum, choosing 
a normalization constant which minimizes x^- The net yield was found to not differ 
beyond statistics using Mi values in the range 1.4 < Mi < 1.8 GeV/c^. Lower values of Mi 
result in over-subtraction and higher values suffer from statistical uncertainties. Note that 
the true pair acceptance for decay of a nuclear resonance with M > 1.4 GeV/c^ is vanishing. 
A value of Mi = 1.4GeV/c^ was selected for the analysis. The solid and dashed histograms 
in Fig. ^ show the results for the true pair invariant mass (M^^^) distribution and the 
normalized combinatorial background distribution. 

The robustness of the analysis procedure was tested via Monte Carlo analysis. A "neg- 
ative test" data sample of uncorrelated pvT^ pairs was generated using the measured single 
particle distributions. Analysis of these data shows no A signal. Additionally, a "posi- 
tive test" was performed on pvT"*" pairs from RQMD-generated events to verify that the A 
resonance is observed correctly there. 

The Miny spectrum of p7r+ pairs after background subtraction is shown in the bottom 
panel of Fig. ^. The total yield into the E814 acceptance for 4.01 x 10^ central collisions is 
587 ± 165 A++. We have computed the acceptance of the E814 spectrometer for A"*""*" using 
GEANT. Over the rapidity interval 1.9 < y < 3.1 covered by the experiment the acceptance 
varies between 3 x 10^'^ at y=2 and 1.2 x 10^^ at y=3. Assuming that the shape of the 
A rapidity and pt distribution is close to that measured for protons [0,|l3]) a-^d taking into 
account a track reconstruction efficiency of 73% per particle, the yield corresponds to 1.7 ± 
0.5 A++ per central Si+Pb collision into 1.9 < y < 3.1. Using the A distributions in y and 
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Pt from RQMD gives a very similar result. 

The predicted yield at freeze-out from the RQMD model is in good agreement with 
the measurement. To best approximate our experimental conditions, we have chosen to 
analyze RQMD events in a similar manner as the actual data and extract the predicted A 
yield using a combinatorial mass spectrum and a mixed-event subtraction. The RQMD pre- 
diction is then 1.6 ± 0.3 per event in the rapidity interval 1.9 < y < 3.1, in remarkable 
agreement with our measurement. The same calculation yields a rapidity integrated yield 
of A+^ of 14.7 ±0.9 and finds 35 % of all nucleons in the A resonance at freeze-out. 

Assuming thermal equilibrium, the measured A abundance yields information |]15| on 
the freeze-out temperature of the system, independent of analyses of the slopes of particle 
spectra. To determine the freeze-out temperature T, we calculated the number density of all 
nonstrange baryonic resonances with masses less than 2 GeV/c^ as a function of temperature. 
The calculation closely follows that of [|15[ but takes, in addition, into account the widths 
of all states. The results for population ratios (which are essentially independent of baryon 
chemical potential) are presented in Fig. ^ as a function of T for all included baryons. 
Using the population ratio for A's of 0.36±0.05 determined from our data we extract the 
freeze-out temperature to be T = 138li8 MeV. 

In summary, we have shown that pion spectra from Si-nucleus collisions at AGS energy 
exhibit a significant enhancement at low pt. A simple model incorporating pions from A 
decay accounts quite accurately for the observed shape if A/nucleon ratios in the range 0.36 
± 0.05 are assumed. RQMD calculations are consistent with our result and also reproduce 
accurately the measured pion spectra. Additionally, the results from the analysis of spectral 
shapes are consistent with a direct measurement of the A"*"^ abundance for Si+Pb, which 
yields 1.7 ± 0.5 A++ per collision in the rapidity interval 1.9 < y < 3.1. Our measurement 
thus quantitatively establishes the importance of the A resonance to the dynamics of the 
collision. We find that the same concentration of freeze-out A's simultaneously explains our 
pion enhancement and the directly measured yield of A"*""*". Finally, the results imply that in 
central Si+Pb collisions a fireball is formed with substantial excitation of A baryons which 
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freezes out at T = ISSlH MeV. 
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FIGURES 

FIG. 1. Pion transverse mass spectra for central {cr/ageo = 2%) Si + Pb collisions in different 
rapidity intervals. Starting with y=4.7, the distributions in each successively lower y bin have been 
multiplied by increasing powers of 10. Shown in the inset are the ratios: Data/Boltzmann (points) 
and mt exponential/Boltzmann (dotted line). For more details see text. 

FIG. 2. Pion transverse mass spectra plotted as a ratio of the data to the best fit Boltzmann 
distribution {mt > 300MeV/c^). Also shown are predictions from a model containing direct pions 
and A decay pions in the ratios 'k c^j'K direct = 0.4 and 0.6 (lower and upper curves respectively). 

FIG. 3. Reconstruction of the A++ resonance. Shown in the top panel are the invariant mass 
spectra of p7r+ pairs from the same event (solid histogram) and of p7r+ pairs from mixed events 
(dashed histogram). The bottom panel shows the difference in the true pair and mixed event pair 
spectra in which the A resonance is visible. 

FIG. 4. Thermal occupation probabilities of non-strange baryons with mass m < 2 GeV/c^ as 
function of temperature. For details see text. 
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